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We investigate the eflect of electronic correlations on the transmission phase of quantum coherent 
scatterers, considering quantum dots in the Coulomb blockade regime connected to two single- 
channel leads. We focus on transmission zeros and the associated yr-phase lapses that have been 
observed in interferometric experiments. We numerically explore two types of models for quantum 
dots: (i) lattice models with up to eight sites, and (ii) resonant level models with up to six levels. We 
identify different regimes of parameters where the presence of electronic correlations is responsible 
for the increase or the decrease of the number of transmission zeros vs. electrochemical potential 
on the dot. However, we show that interaction effects cannot reproduce the universal behavior 
of alternating resonances and phase lapses, experimentally observed in many-electron Coulomb 
blockaded dots. Our numerical results strongly suggest that the main experimentally observed 
features are captured by the theory for chaotic ballistic dots of Molina et al., [Phys. Rev. Lett. 
108, 076803 (2012)] incorporating one-particle wave-function correlations but ignoring many-particle 
electronic correlations. 
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I. INTRODUCTION 

Quantum coherent effects in electronic transport such 
as Aharonov-Bohm (AB) conductance oscillations, weak 
localization, and universal conductance fluctuations-'^ 
originate from interferences between partially scattered 
electronic waves. In contrast to classical transport, 
quantum transport is thus fundamentally influenced by 
scattering phases, and the transport properties of elec- 
tronic nanodevices operating at low temperatures are de- 
termined by complex transmission amplitudes, instead 
of real transmission probabilities. However, while the 
squared modulus of the transmission is the key ingredi ent 
of the Landauer-Biittiker formula for the conductance,!^ 
transmission phases themselves cannot be directly mea- 
sured. 

In their pioneering phase-sensitive experiments, Ya- 
coby et al.'^ measured the conductance oscillations of an 
AB interferometer with a quantum dot (QD) embedded 
in one of its arms. The QD operates in the Coulomb 
blockade (CB) regime and therefore only a single trans- 
verse channel participates in transport. The transmission 
of this channel through the QD is characterized by the 
complex amplitude t = |i|e'". Varying the voltage Vg 



on a nearby plunger gate capacitively coupled to the QD 
allows to add electrons one by one, and the phase of the 
AB conductance oscillations can be monitored as a func- 
tion of Vq . In the two-terminal setup of Ref. [3| , sketched 
in Fig. [1] the conductance reads 

ffAB(Vb,0)=5i°j^(^G) 

+ E-9AB(^G)cos(^27rp^-|-;3,(T/G)), 

where is the flux through the AB ring and — hc/e 
is the flux quantum. The conductance exhibits AB oscil- 
lations with Vc-dependent characteristic phases j3p. An 
Onsager reciprocity relation dictates that g is an even 
function of in a two-terminal setup.'^ Thus, /3p(Vg) = 
0, TT are the only two possible values.'S' In the experiment 
of Ref. g] , /3i was monitored, and abrupt jumps between 
those two values were observed at values of Vg corre- 
sponding to CB resonances, i.e. where an electron is 
added on the QD. Assuming that /3i is directly related 
to the transmission phase a, such jumps could be ex- 
plained by Friedel's sum rule. More puzzling, however, 
were the additional, equally abrupt jumps of tt systemat- 
ically observed in the CB conductance valleys in between 
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FIG. 1; Aharonov-Bohm interferometer, threaded by a flux 
cj>, with a quantum dot embedded in its upper arm. The 
dot is capacitively coupled to a plunger gate potential Vg, 
which tunes its electrochemical potential, and thus changes 
its electronic occupancy. The unsymmetrical shape of the 
dot reflects the presumed lack of geometric symmetry in ex- 
perimental dots. The dark segments at the entrances of the 
dot represent tunnel barriers. The dashed lines on the arms of 
the interferometer stand for a number of possible additional 
leads. 



each and every two consecutive CB resonances. 

Two fundamental questions have been raised at that 
point. First, what is the connection between the conduc- 
tance phases /3p and the transmission phase a? In other 
words, under which conditions is it possible to extract the 
transmission phase a from the experimentally measur- 
able phases /3p? Second, what is the physical mechanism 
responsible for the in-phase behavior, i.e. the systematic 
TT-jumps observed in between CB resonances? 

It was understood earlj^^ that the two-terminal setup 
had to be abandoned to probe the transmission phase 
a. Opening the system to more terminals lifts the reci- 
procity constraints and allows for a one-to-one corre- 
spondence between a and This was experimentally 
achieved by Schuster et alW who opened the arms of the 
interferometer to additional grounded terminals - this is 
sketched by the dashed lines in Fig. [l] Working with such 
a "leaky" interferometer suppresses processes with mul- 
tiple windings around the ring, so that only /3i can be ex- 
tracted. An appropriate tuning of the opening of the ring 
arms in this multi term inal setup allows for the identifica- 
tion of /3i with and Ref. [5] obtained the expected 
Breit-Wigner behavior of the phase, with a smooth in- 
crease of TT every time a CB resonance is crossed. How- 
ever, the second puzzle of Ref. 4; persisted, as the sys- 
tematic phase lapses of tt in-between any pair of consec- 
utive CB resonances also appeared in the multiterminal 
setup. 

Refs. [H and [5] both work with hundreds of electrons 
on the QD. In a more recent experiment, Avinum-Kalish 
et al^ investigated small QD's with zero to few tens of 
electrons. Their key observation is that, as the number 
of electrons on the QD is reduced from 20 down to 0, 
Pi undergoes a crossover from the universal phase be- 



havior regime, with regularly alternating tt jumps at and 
in-between CB resonances, to a mesoscopic regime where 
phase lapses in-between CB resonances occur in a random 
fashion, so that the in-phase behavior of the transmission 
at CB resonances gets lost. These experiments provide 
an important hint towards the resolution of the puzzle: 
candidate theories have to be able to explain the univer- 
sal to mesoscopic crossover as the QD is depopulated. 

This seminal series of works motivated further exper- 
iments. The role of the magnetic field was explored by 
Sigrist et al. in a AB ring with one QD embedded in each 
of its arms.li^ The phase of a QD in the Kondo regime 
was measured by Ji et al^^^^^ Highly controlled experi- 
ments coupled with detailed theoretical modeling for an 
AB device without QD also found phase lapses for some 
parameter values due to scattering and reflections in the 
arms of the ring.li^ 

The puzzles posed by the experimental data attracted 
a sustained theoretical interest. Refs. |7|8|14|15] estab- 
lished that, under not too restrictive constraints, a can 
be extracted from /3i in multi-terminal geometries. As- 
suming therefore that /3i = a, Refs. p!MT^ enounced 
the simple rules that in non- (or weakly) interacting sys- 
tems, TT phase jumps occur under the two following cir- 
cumstances: (i) the electrochemical potential crosses an 
eigenmode of the scatterer (so that one electron is added 
to the QD); (ii) the transmission t vanishes. For non- 
(or weakly) interacting systems, the universal regime of 
transmission phase thus implies that there is one trans- 
mission zero in-between any two consecutive CB reso- 
nances. Numerical simulations on a non-interacting dis- 
ordered diffusive lattice modeF^ however showed that 
transmission zeros occur in-between consecutive reso- 
nances with probability close to 1/2. From this result 
it is often concluded that noninteracting theories are un- 
able to explain the experimental data. 

A mechanism based on level occupation switching ini- 
tially considered in Refs. j20l21j has been repeatedly used 
within non-interacting^^ and interacting models .^^SHsSl A 
Fano-like scenario is assumed to stem from a given 
(broadened) QD level which is much more strongly cou- 
pled to the leads than all nearby (narrow) levels. Then, 
while the position of CB resonance peaks is determined 
by the energy of the narrow levels, the conductance is 
dominated by the transmission through the broadened 
level. As Vq is varied, the narrow levels are successively 
populated right after the CB resonance, at which point 
the broadened level is abruptly depopulated. It was ini- 
tially argued that level occupation switching arises from 
specific spatial structures of the QDj ^° l ^^ lbut later works 
went further and suggested that it generically follows 
from electronic correlations.'^SHU! 

More recently, lattice models for interacting fermions 
were numerically investigated, and an interpretation was 
proposed in which the electronic correlations induce the 
mode switching mechanism. Karrasch et al. investi- 
gated few-level, strongly interacting systems where the 
in-phase behavior is obtained when the single-particle 
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FIG. 3: (Color online) Lattice model for a QD with AI = 8 
sites (red) connected to two external one-dimensional leads 
(black). 



FIG. 2: (Color online) Sketch of the topology of the Resonant 
Level Model. 



level spacing on the QD becomes smaller t han t he level 
broadening due to the coupling to the leads .1^^!^ Varying 
these parameters (and the interaction strength) allowed 
to drive the transition between the universal and meso- 
scopic regime within a given dot of fixed (small) size and 
number of electrons, unlike the experimental case where 
the transition is obtained by changing the electron num- 
ber and thus the filling of the dot. The mode-switching 
mechanism becomes relevant when the electronic popu- 
lation is large enough to justify a mean-field treatment 
of interactions. This conclusion is somehow at odds with 
the numerical results of Ref. [Hj, given that a mean- 
field approximation essentially delivers a single-particle 
theory. Furthermore, the reported behavior of the con- 
ductance (?(Vg) strongly departed from the usual CB 
sequence of alternating conductance peaks and valleys 
found in Rcfs. [4 6 9 lOj, and no arguments were made 
about the generality of the models and parameters con- 
sidered. 

Bergfield et al. considered strongly correlated models 
of molecules, where the universal behavior of the phase 
cannot be reached, unless spatial symmetries are imposed 
on the molecule itself and the molecule- lead couplings.'^ 
This latter result indicates that the universal regime 
requires irregular single-particle spectra, very different 
from regular molecular orbital spectra. Gurvitz proposed 
that the phase behavior in the transmission through a 
quantum dot results from the formation of a Wigner 
molecule.l^ 

In another line of work, a very simple solution to the 
puzzle posed by Refs. [4I6I9| was recently proposed.'^S 
The approach is based on the Constant Interaction Model 
(CIM), which treats QD in the CB regime as noninter- 
acting, up to a constant charging energy term. The sole 
assumption that wave-functions have quantum chaotic 
spatial correlations'^* is able to reproduce the two main 
experimental observations of (i) long, universal sequences 
of in-phase resonances, and (ii) a crossover to a meso- 
scopic regime for not too large number of electrons on 
the QD. The probability of deviating from the universal 
behavior can be obtained as a function of the electron 
filling. Given the success of the G IM in describing dif- 
ferent features of GB physics it was expectable that 
the statistical behavior of the transmission phase were 



also within its reach. 

To bring the experimental-theoretical controversy to 
closure, it is nevertheless important to develop and in- 
vestigate more realistic models incorporating electron- 
electron interactions beyond the GIM (in particular in- 
cluding strong correlations), and check if they also can 
reproduce the experimental observations, for instance via 
the level occupation switching mechanism. This is one of 
the main goals of this work. In particular, we are in- 
terested in knowing if the occurrence of transmission ze- 
ros obtained with particular models of interacting elec- 
trons can account for a genuine universal regime with 
in-phase behavior over long sequences of resonances. For 
that purpose we use two different, generic models of QD: 
the resonant level model (RM), describing a finite num- 
ber of single-particle states with repulsively interacting 
electrons, individually connected to external Fermi liquid 
leads (see Fig. [2]), and lattice models (LM) with nearest- 
neighbor electronic repulsion, connected to two Fermi 
liquid leads via two sites (see Fig. [3|. Our investiga- 
tions are restricted to the experimentally relevant regime 
of zero temperature linear response and we do not ad- 
dress the extreme cases of models with only one or two 

resonan t levels nor situations in which Kondo physics is 
relevant.Eimnol 

We use the DMRG-based embedding method to ex- 
tract the transmission properties of the systerEp3l22l and 
its recent extension to calculate the transmission phase.l^ 
We review such an approach in Appendix [A] and illus- 
trate its power for a one-dimensional interacting quan- 
tum wire, including a numerical verification of the Friedel 
sum rule in Appendix [B} 

The manuscript is organized as follows. In Sec. [ll] we 
introduce the one-particle version of the lattice and res- 
onant models that we use throughout the paper. Follow- 
ing the standard litterature, we relate the transmission 
zeros to a matrix resolvent and to the distribution of the 
coupling widths, which allows to identify the different 
situations in which a transmission zero can appear. Al- 
though such an analysis is not directly applicable when 
correlations are present, it can be extended to the many- 
particle case whenever an effective single-particle theory 
can be constructed. Furthermore, we will present numer- 
ical evidences that the behavior observed in one-particle 
models often also applies to strongly correlated cases. 
The many-particle RM is defined and studied in Sec. 
Ill Results in the limit of large resonance widths are 
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shown, where new zeros appear but no in-phase behavior 
of neighboring resonances is achieved. The next Sec. IV 
is devoted to many-particle lattice models, where exam- 
ples are shown for different sizes and ratios between the 
one-particle level spacing and the level width. No sig- 
nificant trend towards universal behavior emerges upon 
inclusion of correlation effects. Finally, some conclusions 
and perspectives are discussed in Sec. IVl 



II. TRANSMISSION ZEROS IN ONE-PARTICLE 
MODELS 

Lattice models of a one-particle QD connected to one- 
dimensional leads were first used by Levy-Yeyati and 
Biittikci^^ for the numerical calculation of the transmis- 
sion phase. The one-particle resonant level model, where 
the eigenstates of the QD are linked to one-dimensional 
leads through hopping amplitudes, was solved in Ref. [7 
and the necessary conditions for the appearance of trans- 
mission zeros were determined. These two models have 
been extremely useful for discussing the behavior of the 
transmission phase in the non-interacting case, as well 
as when electron-electron interactions are treated within 
the CIM or at the mean-field level. In addition, both 
models can be further generalized to fully account for in- 
teractions by adding a Coulomb term coupling the non- 
interacting basis states. This is the route that we follow 



ti) can be different from the one in the leads, and we use 



in Sections III and IV and thus we start by presenting 
the main concepts concerning one-particle models that 
will be later generalized to describe the interacting case. 
Moreover, our discussion of one-particle models will help 
to systematize and classify the various approaches previ- 
ously proposed to analyze the existence of transmission 
zeros. 

An arbitrarily-shaped, non-interacting QD with M 
sites, connected to one-dimensional leads through its I*^* 
and A/"^ sites, as sketched in Fig. [sj is generically de- 
scribed by the Hamiltonian 

H = Hu + Hg + Hc + Hi^. (2) 
The left and right lead Hamiltonians are given by 



i— — oo i—M-\-l 



(3) 

Cj (cl) indicate the standard operators for annihilation 
(creation) of a spinless fermion on site i and h.c. stands 
for the Hermitean conjugate. We have chosen throughout 
this work the hopping amplitude in the leads to be the 
unit of energy. 

Within the lattice model, the QD Hamiltonian reads 

M 



where (ij) denotes a pair of nearest neighbor lattice sites 
and Hi = cjcj. The hopping amplitude within the dot 



such a freedom in Sec. IV In the absence of magnetic field 
the hopping amplitudes can be chosen to be real without 
loss of generality. Disorder is modeled by random on-site 
energies e^, taken from a uniform distribution of width 
W. 

The capacitive coupling to a nearby gate is propor- 
tional to the total number of electrons in the QD, 
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M 

'Vg ^ Tli 

i=l 



(5) 



where capacitances are included into the definition of Vq- 
Since only the sites i — \ and i — M are connected 
to the leads, the coupling Hamiltonian between the leads 
and the QD is 



TjLM 



-tc (4ci +cIco) -tc[ 



(6) 

The left and right hopping amplitudes connecting the 
leads to the QD could in principle be different from that 
of the leads in order to achieve the regime of weak cou- 
pling. However, we will restrict ourselves to ~ 1 
throughout this work. The choice of one-dimensional 
leads is justified since in the tunneling regime charac- 
terizing CB physics, only a single transverse mode per 
lead is relevant.*^ 

Having defined our one-particle lattice model, we next 
follow Ref. [TS] and discuss the conditions under which 
a transmission zero occurs between two resonances. Our 
starting point is the retarded Green's function of the QD, 
which can be written as 



G(e) = (e - i/D - 



(7) 



with the self-energy S(e) arising from the coupling to 
the leads. The transmission amplitude is related to the 
Green's function via the Fisher-Lee relatiorPil 



t{e)^th{viv,y/^GiM{e), 



(8) 



where ui(r) is the velocity in the left (right) lead for elec- 
trons with energy e. When leads are connected to a single 
QD site, the condition for having t = can be written 



CiM{e-HD + S(e)}-0 



(9) 



where CiAf {A} stands for the cofactor of the (1, M) ma- 
trix element of A. Since Sn and Smm are the only non- 
zero elements of the matrix E, the condition expressed 
in Eq. ^ is entirely determined by the properties of 
the isolated QD and is thus independent of the coupling 
strength to the leads. This important observation allows 
us to locate the transmission zeros from those of the ma- 
trix element 



(10) 
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of the resolvent for the isolated QD, where Cm and 
fpm{i){ipm{M)) are the m"^ QD's eigenenergy and eigen- 
function evaluated at site respectively. 

The structure of ^iM(e) is characteristic of physical 
situations where resonances are coupled to a continuum, 
and allows to determine the existence of zeros according 
to the residues of the poles. As we will see, the reso- 
nant model, that we present below, leads to equivalent 
conditions. 

The one-particle resonant level model is obtained after 
a basis transformation of the QD's degrees of freedom 
from the site basis to the QD's eigenbasis. The Hamilto- 
nian in Eq. ^ retains the same structure, with however 
new QD and coupling terms. 



ttRM 



M 



(11) 



m— 1 

M 



= E Undlco+YmdlcM+i + h.c.) . (12) 



The new fermionic operators dm — X^n ^ni('^)''« 
tained from the old ones via a unitary transformation 



with the eigenfunctions ■i/'m of , and 



S d 

In this way the levels em are understood as eigenval- 
ues of an isolated QD. Alternatively they can be inter- 
preted as the on-site energies of a tight-binding model 
with the topology of Fig. [2] and hopping amplitudes given 
by the partial- width amplitudes Jm. The total widths 

aTeTm = \ll\^ + \Ym?. 

An exact solution for the transmission amplitude in 
this model was presented in Ref. [7J, which obtained 



t{e) = 



2i /ir sin k 



(./n + e-'=)(/„ + e-^-)-|/irP' 
where k is the wave-vector in the lead and 



im ira 



(13) 



(14) 



for X, y = 1, r. With one-dimensional leads the partial 
amplitudes are proportional to the values of the resonant 
wave-functions at the extreme points. Therefore, /ir(e) 
is simply proportional to the function J'iM(e) of Eq. ( 10 ) 



and to the /^- matrix element of the corresponding scat- 
tering problem.'^SMl 

When there are only small variations in the values of 
the wave-functions on the sites connecting to the leads, 
the behavior of J-im{(-) away from e„ is dictated by the 
two surrounding singularities. On the other hand, large 
fluctuations of ';/'m(l)V'm(-M) among different m might 
lead to values of J^ia/ (e) determined by far-away reso- 
nances. These two possible situations will be respectively 
referred to as Restricted Off-Resonance (ROR) behavior 
and Unrestricted Off-Resonance (UOR) behavior. This 
distinction plays a key ro le within the analysis of trans- 
mission phases .fiSHIIlMll] Sj^Qg large wave-function fluc- 
tuations in generic systems are rare, we will find that the 



ROR is the most commonly encountered scenario. In- 
terestingly, the above classification is not only relevant 
for the one-particle models, but also for many-particle 
models (Sees. Ill and IV I. 

The existence of a transmission zero between the m 
and the (m -f I)''* resonances depends on the sign oP^EIl 



th 



V'm(l)V^™(M)V'™+l(l)V^m+l(Af) . (15) 



In the UOR case, when £>,„ > (equal parity of the res- 
onances) there is one zero between the resonances while 
for Dm < (opposite parity) there is no zero. This sign 
rule has been at the basis of several studies of the trans- 
mission phase. In the UOR case we have that for Dm > 
there is an odd number of transmission zeros between the 
m^^ and (to + 1)*** resonances, yielding an accumulated 
phase of tt between the two resonances. For Dm < 
there is no transmission zero or there is an even number 
of zeros, resulting in zero total phase shift between the 
two resonances. In the interferometric experiments on 
QDs operating in the CB regime it is extremely difficult 
to follow the transmission phase across the conductance 
valleys, and only the total phase between resonances is 
relevant. Therefore, the sign rule is also useful in the 
UOR case. At this point it is important to remark that 
the occupation switching mechanism is based on a large 
fluctuation of the partial width leading to the UOR be- 
havior. As we have seen, the appearance of transmission 
zeros when Dm < is indeed possible in the UOR case, 
but they are bound to come in pairs, without altering the 
in-phase relationship of adjacent resonances. 

We next discuss the ROR-UOR competition and the 
sign rule in the one-particle RM. The experience gained 
in this simple one-particle situation will prove very useful 
for understanding the many-particle results in the next 
sections. Compared to the LM, the RM has the advan- 
tage that partial-width coupling amplitudes -f}^ can be 
tuned at will. This flexibili ty has been a key ingredient 
in several theoretical works.'^^J^ ^'^ ' ^^ ' ^^ ' ^" ' 

We first show in Fig. [4] results for a system with four 
equidistant resonances and coupling parameters Jm — 
7j„ = 0.2, for all to = 1,2,3,4, with the exception of 
7| = —0.2. The values of the couplings all have the 
same magnitude, thus leading to ROR behavior. The 
upper panel of Fig. |4] shows the resulting function /h- 
when varying the gate voltage Vg (which is equivalent 
to varying the energy e of the incoming particles from 
the lead) in arbitrary units. The middle panel of Fig. |4] 
shows the transmission coefficient |tp and the phase a as 
a function of Vq. We can see that the presence and the 
absence of zeros between resonances is clearly correlated 
with the sign of Dm- Transmission resonances and ze- 
ros correspond, respectively, to singularities and zeros of 
/if. The trajectory of the transmission amplitude in the 
complex plane as a function of Vg is shown in the bottom 
panel. The trajectories crossing the origin of the complex 
t-plane occupy only one half-plane for the first three res- 
onant peaks, and are characteristic of adjacent resonant 
peaks with in-phase behavior. Even in this simple ex- 
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FIG. 4: (Color online) Numerically obtained data for a 
four-resonance one-particle resonant level model with £,„ = 
m - 2.5, 7™ = Yr,^ = 0.2 for m = 1, 2, 3 and 7! = -7! = 0.2, 
exhibiting ROR behavior in all intervals. Top panel: /ir (de- 
fined in Eq. 14 1. Middle panel: Transmission coefficient \t\^ 



(colored data points) and phase a (solid black line) as a func- 
tion of Vq. Bottom panel: trajectories of the transmission 
amplitude in the complex plane when Vg is varied as on the 
other panels. The color rainbow scale gives the correspon- 
dence between data in different panels. 



ample there is an anomalous extra transmission zero at 
Vg ~ 2.3 far away from the area of resonance peaks. 
This is a finite-size effect, arising from zeros of /xy that 
lie outside the interval where the resonances are concen- 
trated. Such a behavior is irrelevant for CB experiments 
with many resonances, and we only note that, accord- 
ing to our nomenclature, it is an UOR behavior, as the 
value of the transmission far from resonance is dominated 
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FIG. 5: (Color online) Numerically obtained data for a four- 
resonance one-particle resonant level model with em ~ m — 
2.5, 7I = 7i = -0.6 , 72 = 73 = 7l = -0.2, 7^ = -0.05, 
73 — 0.05, and 71 = 0.2, exhibiting UOR behavior between 
the 2"'* and 3'^'^ resonance. The same conventions as in Fig. [4] 
are adopted. 



by the contribution from several resonant tails. Such an 
effect is commonly encountered in finite size numerical 
simulations and it is discussed for instance in Ref. [29] . 

Typical UOR behavior is shown in Fig. [5) where we 
reduced the values of the two central resonance widths 
to make the behavior of /ir between the second and third 
resonances dominated by the external resonances. A loop 
is formed by the trajectory in the complex plane and two 
new zeros appear. These zeros add phase lapses smaller 
than TT in the region between two resonances. The peak 
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appearing around Vq = is related to a small loop close 
to the origin of the i-plane and an accumulated phase 
smaller than it. Therefore, it does not represent a new 
resonance. Instead, it is the result of a broad resonance 
being cut by the extra Fano-like anti-resonances. We will 
see that similar behaviors are also found in the many- 
particle case where the universal behavior is not reached. 

A critical UOR-ROR case arises when the two zeros of 
/ir collapse in a double zero (with a horizontal tangent) 
yielding a kink at the origin of the complex t-plane with 
a phase shift smaller than tt (not shown). This special 
zero is found for specific settings of the parameters and 
if the couplings are perturbed infinitesimally away from 
that singular setting, we find either no zero or two zeros 
in the transmission. 

The inclusion of interactions at the CIM level opens 
a gap of the size of the charging energy in the one- 
particle spectrum between the highest occupied and the 
lowest unoccupied dot levels without affecting the wave- 
functions. The UOR behavior will then be favored since 
many resonances contribute in the electron and hole 
sectors.— However, as explained above, the sign rule dic- 
tating the phase behavior only concerns the two nearest 
resonances. It is therefore important to consider more re- 
fined models including arbitrary couplings and electronic 
correlations. We undertake this task in the forthcoming 
sections and discuss the results in connection with previ- 
ously proposed theories and the concepts introduced in 
this Section. 



III. MANY-PARTICLE RESONANT LEVEL 
MODEL 




-1 I ' ' ' ■ • 1 • 1 

-1 -0.5 0.5 1 



Ret 

FIG. 6: (Color online) Upper panel: transmission coeffi- 
cient (colored) and phase (black) for non-interacting {U = 0, 
dashed) and interacting (U = 2, solid) resonance level models 
with three levels and parameters ei = —0.5, £2 ~ 0.5, ea — 1.5, 
7l = 1.0, = 1-0,72 = -0.5,7^ = 0.5, 73 = 0.4,7^ = 0.7 
Lower panel: trajectory of the transmission in the complex 
plane following the convention of the last three figures. 



The one-particle resonant level model is particularly 
useful to investigate the possibility of UOR behavior, 
because partial widths can be tuned at will. Since in- 
teractions at the CIM level do not favor the universal 
behavior of the phases, it is natural to ask whether cor- 
relation effects beyond mean field could be responsible 
for the experimentally observed universal behavior. This 
question has been analyzed by Karrasch et al.,^^ via 
numerical investigations of many-particle resonant level 
models with up to M = 4 levels and arbitrary level- 
lead couplings. They identified several regimes, deter- 
mined by three energy scales: the mean level-width F, 
the mean one-particle level spacing A in the QD, and 
the strength of the electron-electron interactions U. In 
the case F < A, interactions were not found to affect 
the behavior of transmission phases. In contrast, when 
[/ ~ F > A, sufficiently strong correlations appeared to 
favor the appearance of additional transmission zeros not 
predicted by the sign rule. This result was interpreted 
as the signature of Fano-like antiresonances between a 
renormalized wide resonance with several narrow single- 
particle levels, thus providing some justification for the 
level occupation switching mechanism. While the phase 
lapses reported in Ref. [291 seem consistent with the uni- 



versal regime, the gate-voltage dependence of the trans- 
mission is totally different from the experimentally ob- 
served standard CB resonances. This is due to an incom- 
plete filling of the resonances in the theoretical model re- 
sulting in an associated phase accumulation smaller than 
TT, as we demonstrate below. 

To better understand correlation effects on the trans- 
mission phase we reproduced and extended the results of 
Ref. [29], working with up to M = 6 levels in the QD. 
For M = 2 and 4 our numerical calculations using the 
embedding technique (with Density Matrix Renormal- 
ization Group calculations) agree with those obtained in 
Ref. [29 with the Numerical Renormalization Group and 
the Functional Renormalization Group algorithms. The 
QD Hamiltonian of the many-particle resonant model is 



that of Eq. 11 plus an interaction term 



rrRM 



u 



(16) 



In Fig. |6] we show that interactions can generate trans- 
mission zeros in dots with F > A already for M = 3 
states. The fiuctuations of the level-widths are small 
enough that the non-interacting case (dashed lines) ex- 
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hibits pure ROR behavior. According to the sign rule, 
and since there are alternating parities between consecu- 
tive resonances, there is no transmission zero. The broad 
peak (around Vg = —1) of the non-interacting case (top 
panel, dashed lines) is the result of the overlap of two 
nearby resonances. In the interacting case with U — 2 
(top panel, solid lines), three peaks (around Vq = 0, 3, 
and 5) can be seen. Two zeros are clearly identified from 
the phase lapses of a and from the trajectories of t in the 
complex plane (bottom panel). The small peak between 
the two zeros corresponds to a small loop in the complex 
t-plane, indicating the emergence of UOR behavior in- 
duced by interactions. The increase of a through each of 
the peaks is smaller than tt, which would be the expected 
value if one electron were added to the QD. 

Our example is consistent with the mechanism put for- 
ward by Karrasch et al. |29j , where the transmission am- 
plitude in the vicinity of a new zero is not determined by 
the two nearest resonances but by an anomalously broad- 
ened level. However, the transmission phase change from 
one of the new zeros to the next one docs not correspond 
to the full addition of an electron on the QD. As a direct 
consequence, the Vc-dependence of the conductance does 
not exhibit the characteristic CB peak structure. We ob- 
served similar behavior in all the realizations where extra 
zeros appeared with interactions. The analysis of these 
simple cases illustrates the usefulness of the discussion 
presented in the previous section. Even if the function 
fh{s) of Eq. ( 14) does not have meaning in a many-body 
situation, the concepts of ROR and UOR behavior can 
be addressed by varying Vq and studying the resulting 
trajectory in the complex t-plane. 

The advantage of our DMRG-based embedding 
method is that it allows to increase the QD size to larger 
values than previously studied, still keeping a good pre- 
cision. We next extend the above analysis to larger sys- 
tems, to investigate the possible crossover between the 
mesoscopic and the universal regime of the transmission 
phase. In all examples to be discussed, we randomly 
chose the widths 7™ , while tuning the values of the res- 
onance energies Cm in order to achieve different behavior 
in particular examples. 

We show in Fig. [7] data for a QD with M = 6 resonant 
levels, where the average energy spacing is much smaller 
than the average level width A ^ F. The fluctuations of 
the coupling amplitudes are not strong enough to induce 
UOR behavior in the non-interacting case. The sign rule 
in the ROR case dictates the absence of transmission 
zeros between the resonances (upper panel) , and there is 
one zero outside the region of the resonances. For U = 2 
(central panel) there is one transmission zero outside the 
region of the resonances and four zeros in the region of 
the resonances. In the same way as the previous example 
of a smaller QD, the phase evolution does not exhibit in- 
phase behavior from one resonance peak to the next, due 
to the non-integer filling of the dot at each resonance. 
The trajectory of the transmission in the complex t-plane 
shows loops typical of UOR behavior. For instance, the 
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FIG. 7: (Color online) Transmission coefficient (color) and 
transmission phase (black) of a RM with six levels for (7 = 
(upper panel, dashed) and U = 2 (center panel, solid lines 
and data points). The model parameters are ei = —0.113, 
£2 = -0.071, £3 = -0.059, £4 = 0.004, £5 = 0.016, ea = 0.082, 
7I = 0.3883, Yi = 0.1775, 72 = -0.2106, 72 = 0.3292, 
7;^ = 0.2448, 73 = 0.1214, -ji = 0.2581, 7J = -0.1936, 
75 = 0.3155, 75 = 0.2332, 7^ = 0.4253, 76 = -0.1977. 
Bottom panel: trajectories of the transmission amplitudes in 
the complex plane as a function of Vg- The insets show Vg- 
intervals where transmission zeros occur together with over- 
lapping resonances. 



two resonances close to Vg = are associated with a 
phase increase larger than tt, and are "cut" by a zero at 
Vg — 1- The insets in the lower panel show loops between 
two zeros that are characteristic of conductance peaks 
that do not represent a resonance with the corresponding 
integer filling of the dot. The phase evolution in the loop 
is then smaller than tt, and the angle of the crossing of 
the trajectories at the origin directly gives the missing 
dot filling. 

We next explore the emergence of UOR behavior by 
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FIG. 8: (Color online) Transmission coefficient (color) and 
transmission phase (black) of a RM with six levels for U = 
(upper panel, dashed) and (7 = 2 (central panel, solid). 
The coupling amplitudes are the same as in Fig. [7| while 
the one-particle energies are ei = —1.7497, £2 ~ —1.0161, 
£3 = -0.6483, £4 = -0.4684, es = 0.0428, and ee = 0.16212. 
Bottom panel: Trajectory of the transmission amplitude in 
the complex plane as a function of Vg- 



taking different ratios of F/A within the various reso- 
nances of a given sample, considering the influence of the 
interaction strength U. We present results for a QD with 
M = 6 resonant levels, where the energy spacing between 
the upper two levels is much smaller than all other spac- 
ings for different interaction strengths: [7 = and 2 (Fig. 
[S]) and U — 4 (Fig.joj). For U ~ (dashed Hues) we are in 
the regime of overlapping resonances. This is clearly indi- 
cated by the transmission probability in the top panel of 
Fig. [8j There is an anomalous transmission zero outside 
the resonance area, arising from the same mechanism as 
discussed in Sec. [ll] The sample shows ROR behavior at 
[/ = 0, and in agreement with the sign rule, there is no 
transmission zero. When the interactions are turned on, 
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FIG. 9: (Color online) Results for a resonance level model 
with six levels and the same parameters as in Fig.jS] at the in- 
teraction strength 11 = 4:. Top panel: Transmission coefficient 
(colored line/data points) and phase (black line). Bottom 
panel: Trajectory of the transmission in the complex plane as 
a function of Vg- 



the charging energy leads to the separation of resonances, 
which evolve into CB peaks. Resonances no longer over- 
lap, except the last two, whose energy separation in the 
non-interacting case is particularly small. The change in 
the behavior of the trajectory of the transmission ampli- 
tude in the complex plane induced by the effect of the 
interactions is remarkable. Two new zeros appear rather 
close to each other between the second and the third res- 
onance thus exhibiting UOR behavior with a vanishing 
accumulated phase shift (see insets of the lower panels) . 
They are located in the region where the phase evolution 
is driven by the beginning of the filling of the third level. 
No new zeros appear between the most overlapping last 
two resonances. As the interaction is increased further 
the qualitative behavior and the number of zeros is not 
changed. We show in Fig. |9] the results for stronger in- 
teraction strength, U — 4, where the main difference is 
that CB peaks become narrower and more distant. The 
two zeros between the second and the third CB peak also 
become more separated by the effect of the increasing in- 
teraction. We checked for J7 = 6 and 8 (not shown) that, 
except for this trivial effect, the trajectories of t in the 
complex plane almost do not change with U. 

In a last example of the resonant model, we kept the 
same couplings as in the previous examples, while sig- 
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FIG. 10: (Color online) Top panel: Transmission coefficient 
(colored line) and phase (black line) for a resonance level 
model with six levels and U = 2. The coupling amplitudes 
are the same as in Fig. [7j while the one-particle energies are 
= -0.44976, ez = -0.2161, £3 = -0.04683, £4 = 0.0042, 
£5 = 0.01621, and ea = 0.018. Bottom panel: Trajectory of 
the transmission in the complex plane as a function of Vg- 



nificantly reducing the level spacings to enter the regime 
A < r of Ref. [15]. We see in Fig. fTol that the number 



of zeros does not change with respect to the interacting 
cases of Figs. [8] and |9] but we can obtain very small loops 
in the i-plane, resulting from the extreme UOR behavior. 

The previously discussed examples show that the ex- 
trapolation towards larger systems by increasing F/A is 
delicate, and large variations of this ratio are needed in 
order to generate transmission zeros between each pair 
of consecutive resonances. 

From our numerical results in the many-particle RM 
we conclude that (i) level occupation switching induced 
by interactions appears only in the extreme A ^ F case 
where the zeros can cut through resonances, and (ii) it 
is directly related to UOR behavior where the transmis- 
sion between resonances is not simply given by the inde- 
pendent contributions from the nearest two resonances. 
There is a wide intermediate regime A w F where UOR 
behavior appears only for part of the resonances. The in- 
teraction must satisfy J7 3> A and t/ > F for the effect to 
appear,'22l but once this is the case, increasing U further 
does not induce any qualitative changes in the number of 
transmission zeros and the associated phase lapses. Fi- 
nally, the occurrence of this mechanism should translate 



into a partial occupation of the dot when the special ze- 
ros induced by correlations appear. In this case, the form 
of the conductance peaks significantly departs from the 
standard, symmetric CB peaks. Such a behavior is at 
odds with the experimental observation of well resolved 
CB resonances. 



IV. MANY-PARTICLE LATTICE MODELS 

While the RM is very useful to discuss the various pos- 
sible scenarios that may lead to the appearance of trans- 
mission zeros and the associated phase lapses, the single- 
level energies and half-width amplitudes are independent 
parameters of the model. We therefore return to lattice 
models in this section, where the single-level parameters 
are determined by the geometry of the lattice system. 

The dot Hamiltonian for a lattice model of spinless 
fermions with nearest-neighbor interaction is given by 
that of Eq. Q plus the interaction term 



LM 
DU 



1 



1 



(17) 



Such a nearest-neighbor repulsion can lead to strong 
correlations depending on the value of the interaction 
strength U. Writing the interaction term in the basis of 
the one-particle dot eigenstates, it would take the form 
of (16), but with level-dependent interaction strengths. 



Therefore the many-particle LM is not only less arbitrary 
than its RM counterpart, but also more realistic for de- 
scribing small QDs. The limitation to nearest-neighbor 
interaction considerably simplifies the numerical work, 
but it is not of fundamental nature. 

When U 7^ 0, it is not possible to make a priori general 
statements about the behavior of the zeros. We therefore 
performed numerical calculations using the embedding 
method presented in detail in the App.|A] We considered 
systems ranging from the simplest topology of a diamond 
with M — 4 sites (Fig. 1 of Ref. ,40] ) to dot sizes up to 
M = 8 sites (see Fig.[3|. 

When the single-particle level-spacing is of the order 
of the level- width (A > F), the transmission zeros with 
U typically exhibit the same qualitative behavior as 
in the non- interacting case U = |29|40j . This situation 
typically happens in small systems when <d = 1. For 
instance, in the case with M = 4 and on-site energies 
chosen such that one has a transmission zero at U = 0, 
the effect of the interactions is to separate the peaks, 
reducing their widths without changing the qualitative 
behavior of the scattering phasd^. 

In simple models of small dots without disorder, inter- 
actions do not modify the number of transmission zeros 
from the non-interacting case and they do not induce 
a change from ROR to UOR behavior. By varying the 
on-site energies and taking < 1, we occasionally ob- 
tain the displacement of a transmission zero from the re- 
gion between the resonances to the zone outside the reso- 



nances. Such a behavior is shown in Fig. 11 for the case 
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FIG. 11: Transmission coefficient vs. gate voltage of a lattice 
model for a quantum dot with four sites and tu = 0.25 for 
(7 = (4) in top (bottom) panel. The chosen disorder has 
led to on-site energies ei = 0.220, £2 = 0.360, es = 0.219, 
£4 = —0.037. The logarithmic scale used in the insets allows 
to locate the transmission zeros. 



M = 4, tu = 0.25, U — A and a realization of the on-site 
energies randomly chosen in the interval [— VF/2, W/2] 
with W = 1. In this realization, there is a transmission 
zero between two groups of resonances in the case U = A; 
tD = 1 and also in the non-interacting case U = when 
tD — 0.25. However, this zero disappears in the shown 
case of tu = 0.25 and U = 4, while another zero appears 
outside the resonance area. The disappearance can be 
related to a change in the wave function of the dot that 
evolves into a charge density wave in the interacting case. 

According to Ref. [29], the evolution from the meso- 
scopic to the universal regime can be achieved if A < 
F < [/. To check this idea, we reached the previous con- 
dition within the many-particle LM by going to relatively 
large systems (M = 8) and taking the hopping within the 
dot smaller than in the leads. 

Different sizes (M = 4 to 8) level spacings {tu = 1 
and < 1) and interaction strengths {U = to 8), with 
and without disorder, have been explored in a consider- 
able numerical efforlP^. Below we present some illustra- 
tive examples supporting our general conclusion: inter- 
actions do not typically induce the transition from the 
mesoscopic to the universal regime. 

In general, the main effect of the interaction is to sep- 
arate the resonance peaks and to make them narrower. 
In the more common scenario the number of transmis- 
sion zeros is unchanged from the non-interacting to the 
interacting case. In some cases, the structure of the res- 
onances may substantially change when interactions are 
switched on and the position of the zeros can change 
accordingly. Resonances that are very close and then in- 
distinguishable in the non-interacting case, can often be 



FIG. 12: \t\ vs. Vg of a lattice model for a quantum dot 
with M — 8 sites and U — (A) in the top (bottom) panel. 
tu ~ 0.5 and W — 1. The inset shows the same data on a 
logarithmic scale. 



resolved in the interacting case. A scenario we have ob- 
served in several examples is the disappearance of trans- 
mission zeros beyond a given interaction strength. This 
situation arises because the original UOR zeros occurring 
in the regime F ^ A disappear as interactions change 
the ratio between F and A. Pushing a transmission zero 
outside the interval where the resonances are (like the ex- 
ample shown in Fig. 11) is another common outcome of 



including interactions. All cases analyzed resulted in the 
mesoscopic regime with a random alternation between 
resonances and transmission zeros and significantly less 
transmission zeros than resonances. In no instance did 
we find that interactions can induce a transition from 
ROR to UOR behavior as it was the case in the RM. 
Scaling up the system size or varying tu and U did not 
change these conclusions. 

In Fig. [12] we show an example with AI — 8, U — 2, 
and tu = 0.5 exhibiting the most commonly encountered 
behavior where the number of zeros does not change from 
the non-interacting to the interacting case. Even if the 
structure of the resonances and the position of the zeros 
change dramatically from one case to the other, the total 
number of zeros remains constant. 

We finally illustrate in Fig. JTS] a scenario we often ob- 
served, with less zeros in the interacting case than in 
the non-interacting one in a quantum dot with AI =^ 8, 
U = 2, W = 1, and tu = 1. The double zero in the 
middle of the curve (around Vq — 1) of the U = case 
corresponds to an UOR case that is transformed to a 
normal ROR zero for [/ = 2. As resonant peaks become 
narrower and more separated in the presence of interac- 
tions UOR behavior is not favored. This is very different 
from the situation in the RM, and we do not observe here 
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nances. In the ROR case we have one transmission zero 
or none, depending on the sign of Dm [see Eq. (15l] be- 



FIG. 13: |tp vs. Vg for [/ = (2) in top (bottom) panel, 
of a lattice model for a quantum dot with eight sites, to = 1, 
and W = 1. The inset shows the same data on a logarithmic 
scale. 



the population switching mechanism. 



V. CONCLUSIONS 

In this work we have critically examined the influence 
of electronic correlations towards the observed universal 
behavior of in-phase resonances when a quantum dot op- 
erating in the CB regime is placed in one of the arms 
of an AB interferometer. In order to compute the com- 
plex transmission amplitude of the scattering matrix cor- 
responding to strongly correlated systems we have ex- 
tended the embedding method previously used for con- 
ductance computations. For the simplicity of the numer- 
ical calculations we have used a model of spinless elec- 
trons. This is appropriate since spin effects are only ex- 
pected to be important for extremely small dots, while 
our interest is to achieve the conditions of dots large 
enough to reach the transition from the mesoscopic to 
the universal regime. 

We have analyzed two different models including 
electron-electron interactions: the resonant level model 
and lattice models with nearest-neighbor interactions. 
The non-interacting limits of these models were also dis- 
cussed in order to define different scenarios for the occur- 
rence of transmission zeros and phase lapses. When the 
width of the resonances is smaller than the level spacing, 
the peaks are well resolved in energy. Two cases should 
then be considered: restricted off-resonance (ROR) be- 
havior, in which the two nearest resonances determine 
the character of the transmission amplitude in-between; 
and unrestricted off-resonance (UOR) behavior, where 
the transmission amplitude between two adjacent reso- 
nances is significantly affected by other far-away reso- 



ing positive or negative, respectively. In the UOR case, 
and when the resonances are well separated, the num- 
ber of zeros might be increased with respect to the ROR 
case by a multiple of two, which would not change the 
total phase shift accumulated in the interval between the 
resonances. Therefore, we have a phase shift of tt or 
between the resonances depending on Dm being positive 
or negative, respectively. 

The situation changes when the resonances are 
strongly overlapping and it has no meaning to treat peaks 
and valleys separately in the energy dependence of the 
transmission amplitude. In these cases the transmission 
is typically given by the contributions from many reso- 
nances resulting in an UOR behavior with the possible 
appearance of zeros. However, the phase accumulated in 
the region of overlapping peaks is not an integer multi- 
ple of TT. The transmission zeros cut the phase evolution 
leading to the incomplete filling of the dot and phase 
shifts smaller than tt. This case is clearly not represen- 
tative of the experimental situation, where well-resolved 
peaks alternate with conductance valleys. 

In the case of the resonant model we confirm previ- 
ous result^2ll about the possible increase in the number 
of zeros due to the presence of interactions in the regime 
of very wide resonances, depending on the choice of the 
coupling parameters. We relate this phenomenon to the 
UOR behavior studied in the non-interacting case. How- 
ever, results for the RM do not reproduce key features of 
the experiment Hke the in-phase behavior of consecutive 
resonances. 

We treated lattice models of dots with up to M = 8 
sites. In most of the cases the number of transmission 
zeros was independent of the interaction strength, while 
in a minority of cases we have observed that interactions 
can induce ROR behavior from UOR behavior as res- 
onances become more narrow and isolated, thus reduc- 
ing the number of transmission zeros. Consistently with 
the experimental findings, we observe the mesoscopic be- 
havior for the small size QD that we treat numerically. 
In addition, when decreasing the internal hopping am- 
plitude to achieve smaller energy level separation with 
respect to the level couplings, no tendency towards uni- 
versality was obtained, independently of the value of the 
interaction strength. The exploration of a large parame- 
ter space of level separations, coupling widths, and level 
spacings allowed to approach the conditions of not so 
small QD, where the universality was claimed to arise by 
the effect of electronic correlations. This is not the case. 
Only at large enough QD {kLs ^ 1), the one-particle 
wave-function correlations provoke the emergence of uni- 
versality, but for those relatively large sizes the electronic 
correlations are no longer important. 

In the non-interacting case with chaotic underlying dy- 
namics, the generic distribution of eigenstates and partial 
width&2Si favor the ROR behavior. Taking interactions 
at the CIM level may lead to the UOR behavior, but 
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with well separated resonances. Therefore, in these cases 
the sign rule based on the one-particle wave-functions 
determines the phase behavior. The case of UOR with 
overlapping resonances is achieved by some tuning of the 
system parameters in the non-interacting and CIM cases. 

The main conclusion of our work is that the present 
understanding of the behavior of the transmission phase 
through quantum dots in the presence of strong corre- 
lations cannot by itself explain the observed universal 
behavior for large dots. This result is consistent with 
the observation that the emergence of the universal be- 
havior can be obtained taking into account one-particle 
wave-function correlations but not many-particle elec- 
tronic correlation effects. 
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Appendix A: Embedding method for the 
transmission phase 

The embedding method is a powerful technique to 
calculate the conductance through a strongly correlated 
nanosystenp2it39| ^jth or without disorder. The system 
of interest is connected to a one-dimensional lead that 
closes into itself resulting in a ring which is pierced by 
a magnetic field. The response of ground state proper- 
ties to such a perturbation, like the persistent current or 
the phase sensitivity, allows to infer the conductance of 
the original system. While the one-dimensional setups 
with spinless electrons have been the most commonly 
used models, quasi-one dimensional leads and Hubbard 
chains have been recently considered,'^ and nanosystems 
with non-trivial structure have also been studied.ESl The 
embedding method has been recently generalized in order 
to obtain not only the conductance, but also the trans- 
mission phas^^, and detailed tests of its precision have 
been performed. 

In this appendix we present the main steps of such a 
derivation in order to settle the notation and the basis of 



the numerical method used in Secs. |III| and IV The setup 
of the embedding method is given by a Hamiltonian as 
the one of Eq. ^ , where Hd stands for the Hamiltonian 
of the quantum dot depending on the model considered. 
Hq ^ allows for the application of a gate voltage, while 
the coupling term is given by ([6]) or ([12]). The lead Hamil- 
tonian needs to be modified with respect to ^ in order 
to represent a ring pierced by a flux (p^. It reads 



= - 



L 

E 

i=M+l 



with the boundary condition cq = exp{2TTi(j)(./4'o)cL. In 
the inset of Fig. [M] we show the embedding setup for a 
linear QD (a chain) as used in App. [B] (notwithstanding 
a QD of arbitrary shape like those of Figs. [3] and [2] can 
be treated). Despite the similarity between the embed- 
ding setup and that of the AB interferometer, they are 
very different since the first is a closed system with fixed 
number of particles and the flux 0c is an auxiliary one, 
without physical reality. 

Staying at first within a one-particle model (that is, 
without the interaction terms ( 16 ) or ( 17 )) and with one- 
dimensional leads, the transport through the nanosystem 
is characterized by the scattering matrix 



S = 



r t' 
t r' 



( ze'^ cos Q e"* sin d 
\ e~"' sin 9 ie^*^ cos 6 



(A2) 



We have chosen a generic parametrization of a 2 x 2 uni- 
tary matrix. The transmission amplitude for particles 
coming from the left of the QD and as defined in the 
introduction is given by 

t = |t|e^" ^e^^^-^^sinei. (A3) 

It is related with the Green function by Eq. [s] t' is the 
transmission amplitude for particles impinging from the 
right and r(r') is the reflection amplitude for particles 
coming from the left (right) of the QD. The angle 9 and 
the scattering phase C are restricted to the interval [0, tt), 
while the phases 77 and ^ are defined on [0, 27r). 

When the Hamiltonian of the dot exhibits time- 
reversal symmetry, one has 77 = or 77 = tt. This will 
be our case, since the artificial flux <j>e used to drive the 
persistent current is seen by the ring, but not by the dot. 
When a control parameter is varied (like Vq or e), a jump 
of 77 between its two possible values can only occur when 
the transmission amplitude vanishes (sin 9 = 0), in order 
to preserve the continuity of the scattering matrix. This 
observation is equivalent to that of Sec.|T]about the cross- 
ing of the origin of the complex plane by the transmission 
amplitude being typically associated with a jump of tt of 
its phase a. 

Right-left symmetry would restrict ^ to or tt. Simi- 
larly as in the case of time- reversal symmetry, jumps of 
TT in ^ are only allowed when cos 9 = 0, that is, when the 
reflection amplitude vanishes. However, throughout this 
work, we consider generic QDs with arbitrary ^. 

When considering the phase evolution as a function 
of an external parameter, it is often convenient to work 
with the accumulated phase a^, whose range of definition 
is not restricted to the interval [0, 2tt). 

Embedding the scatterer in a ring of length Ll pierced 
by a dimensionless flux $0 = '^t^4'c/4'q leads to the fol- 
lowing quantization condition for the one-particle states 
of the composed systerrPSl 



cos($p 



V) 



1 



sm( 



cos (fc^L + C) 



(A4) 



(Al) 



Since the scattering phase ( belongs to the interval [0, tt), 
there are two branches for the solutions (in k) of (A4) 
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corresponding to the two possible values of rj (0 and tt). 
On the other hand, the transmission phase a is defined in 
[0,27r), which allows to write (A4) in the more compact 
way 



COs($e) 



1 



sin0 



COS (fcil 



(A5) 



We note Lg the length of the scatterer between the leads, 
i = Ls + Ll the total length of the ring, and the phase 
shift 6a — a — kLg. We express L and k in Anderson 
units (that is, in terms of the lattice spacing). Knowing 
the dispersion relation in the leads £{k), the sum over 
the energetically lowest one-body energies e(fc„) allows 
to obtain the ground-state energy E{N) of the whole 
system containing N particles. For an odd number of 
particles N = 2np + 1, the lowest order terms in a 1/L 



expansion reai 



r, "F 

n=l 



de 
dk 



fo(A:(0)) + O[l/L] 



(A6) 

with the fc-values in a clean ring fci"^ 2Tm/L. The 
ground state energy of a clean ring (with a scatterer hav- 
ing perfect transmission \t\ = 1 and da ~ 0) is given 

by 



£;(o)(7V)=6(0)-f2^6(fc(")). 



(A7) 



71=1 



The embedding method allows to obtain the conductance 
of the scatterer from the flux dependence of E{N), which 
only appears in second order in 1/L [35]. By changing 
the total number of particles, we have access, with the 
help of (A6) and (A7), to the scattering phase shift at 



the Fermi level (fcp = liniL- 



, kn] ) of the lead 



Sa{k-p) 



j.^^^ L ( E{N)-E{N-2)-2e{k^^}) 
L^oo 2 \ de/dfc|^.^^(o) 



(A8) 

In a chain with N ~ L/2 + \ (and L multiple of 4) we 



A8 



are effectively at half filling, kn^ = 7r/2 such that Eq. 
reduces to 



8a{kY) = - lim - [E{L/2 + 1) - E{L/2 - 1)] . (A9) 

L—foo 4 

The previous derivation is based on a single-particle ap- 
proach. As in the case of the embedding method for 
the conductance, the passage to the many-body system 
is justified by the fact that the scattering properties of a 
many-particle scatterer can be represented by an effective 
one-particle scattering matrix. For instance, it has been 
verifiecP^ that the flux-dependence of the many-particle 
ground state is, in the large Ll limit, reproduced by the 
total energy obtained from effective single-particle states 
(as done in the previous derivation). 



In the many-body case, the scattering matrix (A2) is 
understood as an effective one, where each of its enti- 
ties depends on the interaction strength U. When using 



Eqs. or X9 in order to obtain the transmission phase. 



the limiting procedure is implemented by extrapolating 
towards large values of Xl. This procedure is numeri- 
cally demanding, and constitutes the bottleneck of the 
embedding method.'*^ 

The eigenphases of the scattering matrix are, for 
cos cos ^ > 0, ipi = ^ -|- Arcsin(cos0cos^) and Lp2 = 
( + TT — Arcsin(cos 9 cos 5) , and the corresponding Wigner 
time is 



Tie) 



h x ^ dip, 

2^ 



q=l 



de 



de 



(AlO) 



The one-particle density of stat es in t he scattering region 
is related to the Wigner time 



(All) 



where the brackets stand for a spectral average over many 
eigenstates. 

Integrating o ver a n energy interval we recover the 
Friedel sum ruld^^l^ 



AC = vrAiVs 



(A12) 



as a relationship between the number of particles 
added to the scattering region and the corresponding 
change in the scattering phase. 



The fact that the Friedel sum rule in its form (A12) 



applies to the scattering phase C has been emphasized in 
the literatur^iSHlll. The lapses of TT in the transmission 
phase a at the zeros of t are not related with a special 
behavior of the density of states. However, the integra- 
tion leading to (A12) involves a large energy interval (on 



the scale of the level spacing) and therefore the lapses av- 
erage out and the accumulated phase follows the Friedel 
sum rule. 



Aac = ttAA^s 



(A13) 



Such an average behavior has been discussed in Ref. 
|33j where the ambiguity between the lapses of tt and 
— TT was proposed to be lifted by applying a small mag- 
netic field, allowing to avoid the origin in the complex 
i-plane and obtaining well defined phase lapses with a 
phase change close to tt or — tt. The scattering phase can 
then be obtained from the accumulated phase by taking 
[a,{B = 0+) + a,{B = 0-)]/2. 

The existence of phase jumps in the Vo-dependence 
of a can also be obtained from the standard embedding 
method applied to the conductance by locating the zeros 
of the transmission, as we do in Sec. |IV| The quantization 
condition ( A4 ) leads to a phase sensitivity^Si 



AE = E{N, $c = tt) - E{N, $e = 0) 



hvp 
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— — Arccos(sin 9 cos ry) 



and therefore 



\t\ — sin ( AEj^^ COST] — sin 



AE 



A£;(o) 



(A14) 



(A15) 
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where AE'^"^ is the phase sensitivity of a perfectly trans- 
mitting scatterer. Since we are working with time- 
symmetric dots, 77 can only take the values or tt. The 
switches between these two branches may only occur 
when t — 0. Therefore, the sign changes of AE are as- 
sociated with the vr lapses in rj and a. One should also 
notice that there can also be transmission zeros without 
phase lapses in cases when there is a zero of AE with- 
out a sign change. Indeed this possibility can occur for 
particular values of the parameters, as discussed in Sec. 

For a strictly one-dimensi onal s ystem, the sign of AE 
is fixed by Leggett's theorerrPiES] Por an odd number of 
particles we are in the branch of 77 = 0, the transmission 
amplitude t never vanishes, and there cannot be branch- 
switches or phase lapses. For a quasi-one dimensional 
scatterer this is no longer true, and we expect to have 
parameter values where t vanishes and branch-switches 
appear. 



Appendix B: Transmission phase of a 
one-dimensional many-body scatterer 

The applicability of the embedding method to obtain 
the transmission phase of a many-body scatterer can be 
conveniently tested in the one-dimensional case. As dis- 
cussed at the end of [X] one-dimensional systems are con- 
strained to the branch 77 = 0, thus a — and there are 
no transmission zeros. On the other hand, comparing the 
numerical results to the prediction from the Friedel sum 
rule constitutes a valuable test of the method. 

An interacting one-dimensional chain is a particularly 
simple example of a many-particle lattice model where 
the dot Hamiltonian (Eq. |4]) only has M = Lg ordered 
sites. For simplicity, we work in this section in the non- 
disordered case Ci = and we take ~ I. In the ab- 
sence of a gate voltage ((Vq = 0)), we work at half fill- 
ing, E{L/2 + 1) = E{L/2 - 1), for periodic boundary 
conditions in the ring and thus Sa = 0, independent of 
the interaction strength. The number of electrons in the 
scattering region 



(Bl) 



i=l 



is equal to is/2. This is consistent with the findings of 
Ref. ISni , where Fabry- Perot like oscillations of the trans- 
mission through two interacting regions in series were 
studied. 

In contrast, once an additional gate voltage is applied 
(Eq. [5]), particle-hole symmetry is broken and TVs will 
differ from The phase shift is thus expected to be 

nonzero. For such a setup, the embedding method has 
been used to show that Coulomb-blockade like oscilla- 
tions of the conductance as a function of Vg appear in 
the presence of interaction^^ even in the well coupled 
case. The charge in the dot region decreases in steps 




-TT 



-27r 



-Stt 



FIG. 14: (Color online) Accumulated transmission phase shift 
for a many-body scatterer of length Ls = 10 with interaction 
strength U = 2 (solid line), as a function of the gate voltage 
Vg- The squares are the prediction of the Friedel sum rule 
based on the numerically obtained charge densities on the 
Ls interacting sides. The diamonds result when the densities 
on the Ls sites of the scatterer plus that of 10 additional 
sites surrounding the interacting region. The inset shows a 
sketch of the embedding setup for a simple one-dimensional 
interacting scatterer (in red). 



once Vg is increased and depletes the interacting region. 
We extend now these DMRG^'* based calculations to 
compute the properties of rings embedding such a many- 
body scatterer, as well as the transmission phase result- 
ing from (A9). We choose a chain with Ls — 10 sites and 



an interaction strength U = 2. 

Data for ring sizes up to L = 120 were used in the 
extrapolation towards infinite size. The results for the 



transmission phase are presented in Fig. 14 (solid line) 
together with ttA^s — 7rLs/2 (squares). 

Both quantities should be equal according to the 
Friedel sum rule since Sa — a — 7rLs/2 for half filling. 
The results have very similar behavior but small quanti- 
tative differences appear. 

The difference disappears when the density modifica- 
tions outside the Ls sites of the scattering region are 
included in the calculation of A^s- Taking into account 
the density changes in the interacting region plus that 
on 5 additional sites on either side of the scatterer (dia- 
monds), the values of ttA'^s — 7rLs/2 are in quantitative 
agreement with those of 6a. 

From this numerical example we have learned how pre- 
cise the embedding method for scattering phases is, and 
we have checked that in order to comply with the Friedel 
sum rule, all the charge displacement in the neighbor- 
hood of the scatterer has to be accounted for. The effect 
of the charge build-up in the leads of the AB interferom- 
eter was invoked in Ref. [S] as an important ingredient 
in order to address the physics of the experimentally ob- 
served in-phase behavior of consecutive resonances. It 
might be responsible for the phase increase smaller than 
TT at certain resonances which is observed in the experi- 
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mental dataP 

The embedding method is part icula rly efficient when 
dealing with one-dimensional lead^^^I^, but the scatterer 
might have any topology or dimensionality. Upon this 



fact is based our numerical work of Sees. Illll and llVl where 
the quasi-one dimensional systems that allow transmis- 
sion zeros are thoroughly studied. 
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